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ABSTRACT 

We report the results of a pilot multiwavelength survey in the radio continuum 
(X, Ka and Q bands, i.e., from 3.6 cm to 7 mm) carried out with the EVLA in 
order to confirm the presence of very large dust grains in dusty disks and torii 
around the central stars in a small sample of post-AGB objects, as inferred from 
millimeter and sub-millimeter observations. Supporting mm-wave observations 
were also obtained with CARMA towards three of our sources. Our EVLA survey 
has resulted in a robust detection of our most prominent submm emission source, 
the pre-planetary nebula IRAS 22036+5306, in all three bands, and the disk- 
prominent post-AGB object, RVTau, in one band. The observed fluxes are 
consistent with optically-thin free-free emission, and since they are insignificant 
compared to their submm/mm fluxes, we conclude that the latter must come from 
substantial masses of cool, large (mm-sized) grains. We find that the power-law 
emissivity in the cm-to-submm range for the large grains in IRAS22036 is u 13 , 
with (3 = 1 — 1.3. Furthermore, the value of /3 in the 3 to 0.85 mm range for the 
three disk-prominent post-AGB sources ((3 < 0.4) is significantly lower than that 
of IRAS22036, suggesting that the grains in post-AGB objects with circumbinary 
disks are likely larger than those in the dusty waists of pre-planetary nebulae. 

Subject headings: stars: AGB and post-AGB, stars: mass-loss, circumstellar 
matter, accretion disks, binaries: general, radio continuum: stars 
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1. Introduction 

Ordinary stars (i.e., with main- sequence masses 1-8M ) die extraordinary deaths, pro- 
ducing profound effects on their environment before they fade into obscurity as white dwarfs. 
The mass-loss phenomena in these objects play a key role in the Galaxy's chemical and dy- 
namical evolution, and dramatically alter stellar evolution. Post-AGB (pAGB) objects are 
believed to be objects in transition between the AGB and planetary nebula (PN) evolution- 
ary phases and hold the key to some of the most vexing problems in our understanding of 
these very late stages of evolution for ordinary stars. 

We have made considerable progress in understanding the mass-loss processes which 
govern the transformation of the spherical mass-loss envelopes of AGB stars into strikingly 
aspherical PNs (Balick & Frank 2002, Sahai et al. 2011), specifically with the identification 
of fast, collimated outflows as the most likely primary agent for the change in morphology 
(Sahai & Trauger 1998). But the presence of dusty disks (or highly-flattened equatorial 
structures) around pAGB stars remains an enigma: an intriguing feature of these pAGB 
objects is the presence of strong submillimeter excesses, which suggest the presence of very 
large, cold grains. Here we report observations that confirm the presence of large grains 
in two key objects, paving the way for detailed dust-radiative transfer models to accurately 
constrain the grain sizes and theories for the formation and longevity of disks in which grains 
can grow to such large sizes. 

There are two major classes of pAGB objects having equatorially-flattened structures. 
The first class has strong evidence for medium-sized (~50AU) disks and very seldom shows 
prominent nebulosity (hereafter dpAGB, i.e., disk-prominent pAGB objects). The near- 
to-far infrared spectral-energy-distributions (SEDs) of these objects have been modelled as 
arising from disks (de Ruyter et al. 2005, 2006); the disks are believed to be gravitationally 
bound. A large fraction of the central stars of dpAGB objects consists of short-period radial- 
velocity binaries, and many of them show stellar photospheric depletion patterns resulting 
from a poorly understood process in which the circumstellar dust is trapped in a disk, and 
the dust-depleted gas is accreted back onto the star (Waters et al. 1992). Recent near- 
and mid-infrared interferometric imaging with the VLTI of a subsample of these objects has 
provided direct evidence for the disk structures and their sizes (e.g., van Winckel et al. 2008), 
confirming their circumbinary nature. The ISO-SWS (~2.5-40/iin) and Spitzer spectra of 
many of these show amorphous and crystalline silicate features, indicating significant dust 
processing, which is only possible in a stable circumstellar disk, de Ruyter et al. (2006) 
modelled the 0.85 mm fluxes in dpAGB objects as arising from very large (mm-sized) grains; 
detailed 2D disk-modelling by Gielen et al. (2007) has lent further support to this conclusion. 

The second class of pAGB objects consists of pre-planetary nebulae (PPNs), observa- 
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tionally recognized by prominent aspherical nebulosities surrounding a central pAGB star. 
Most bipolar or multipolar PPNs harbor overdense, dusty equatorial waists that are much 
larger than the dpAGB circumbinary disks. In an extensive HST imaging survey of PPNs, Sa- 
hai et al. (2007a) recognize these waists as an important morphological feature in PPNs: for 
many of these, the waists display sharp outer (radial) edges, with radii typically ^ 1000 AU. 
For a few PPNs for which submm data are available, the latter indicate the presence of large 
grains as in the dpAGB objects. For example, in the bipolar PPN IRAS 22036+5306, Sahai 
et al. (2006), find an unresolved (<0'.'85) submm continuum source associated with the cen- 
tral, dusty torus seen in HST images (Sahai et al. 2003). Modelling of the full SED, including 
the submm flux, shows that the latter is produced by a very substantial mass (^0.02M Q ) of 
large (mm-sized), cool (T^50K) grains. Similar results have been inferred for other PPNs, 
e.g., IRAS 18276-1431 and IRAS 19475+3119, from OVRO 2.6mm data (Sanchez Contreras 
et al. 2007, Sahai et al. 2007b), and CRL 2688 from VLA 1.3-3.6 cm data (Jura et al. 2000). 

The connection between dusty waists in PPNe, the disks in dpAGB objects, and inner 
accretion disks is unknown. Yet, these structures appear to be key ingredients in our under- 
standing of the late evolution of most stars - the inner disks are the only ones that are likely 
to be able to contribute to jet production, while the dusty waists are the only ones that 
contribute to the observed morphologies in PPNs. It is believed that accretion disks may 
form from either (a) the shredding of a low mass companion around an AGB primary core 
(Reyes- Ruiz & Lopez, 1999), or (b) Bondi-Hoyle accretion of the primary wind around the 
companion (Morris 1987, Soker and Livio 2004). Smoothed-particle hydrodynamics (SPH) 
codes have demonstrated how an unbound outflow from a mass-losing primary can be concen- 
trated toward the system's equatorial plane into a small, ~(l-few) AU, bound disk around a 
binary companion (Mastrodemos and Morris 1999). The Red Rectangle (HD 44179), which 
is the closest known PPN with a binary central star, provides a connecting link between 
these two classes of pAGB objects - direct evidence for a large (outer radius ~ 300 — 550 
AU), bound disk in this object comes from interferometric CO(2-l) mapping (Bujarrabal et 
al. 2003), and HST imaging reveals a bipolar PPN (Cohen et al. 2004). 

The disk or torus mass can provide an important constraint on the formation process. 
For example, common envelope evolution (Moe & de Marco 2006) would lead to expulsion 
of most of the stellar envelope. Hence, a large value of the disk or torus mass (~0.01M Q in 
dust), as, e.g., in IRAS22036, would support a common envelope origin; wind-accretion 
would lead to much lower values for this mass (as, e.g., found in HD 44179). However, since 
the frequent detection of strong Ha emission (e.g., Maas et al. 2005, Sanchez Contreras et 
al. 2008) in pAGB objects shows that some ionized gas is typically present, the contribution 
of free-free emission from ionized gas to the submm/mm fluxes must be determined before 
reliable masses can be observed. Furthermore, the dust power-law emissivity index in the 
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mm-to-cm wavelength range is poorly known. In order to address both these issues, we 
carried out a pilot multifrequency (43.3, 33.6, & 8.5 GHz) survey of 10 pAGB objects with 
the National Radio Astronomy Observatory'^] EVLA facility and report our results in this 
Letter. 



2. Observations 

We selected 10 post-AGB objects (Tabled]) for EVLA observations on the basis of their 
published mm or submm continuum fluxes (e.g., Gledhill et al. 2001, Giirtler et al. 1996) and 
a new interferometric survey with OVRO (Sanchez Contreras & Sahai 2011) and CARMA 
(reported in this paper). 

The EVLA observations were performed with the array in the D configuration in a series 
of short scheduling blocks under project code AS1021, from March 26, 2010 through June 
21, 2010. The EVLA WIDAR correlator was configured to give two subbands each with 
a bandwidth of 128 MHz and dual polarization. Three frequency bands were observed for 
most sources (though not all, because of time constraints) with center frequencies of: 8396 
& 8524MHz (X band), 33496 & 33624MHz (Ka band), and 43216 & 43344MHz (Q band). 
Observations in the higher frequency bands typically had more time on source in order to 
have similar signal-to-noise ratios in all three bands. For low declination sources, however, 
the increase in the rms noise levels for the two high frequency bands are still apparent. The 
synthesized beam sizes depend upon the receiver band as well as the declination: they range 
from 8"2x2'.'5to 33"3x20"9for X band; 2"2xi'.'8to 9"5 x 1"8 for Ka band; and 1"5 x 1"4 to 
6'.'2 x l'!7 at Q band, depending upon the declination. 

Absolute flux calibration was determined by observing either 3C286 or 3C48, which 
are EVLA flux-density standards. The flux density scale is likely accurate to about 5% for 
X-band, high declination sources, but degrades to 20-30% for Q-band, and especially the 
low-declination sources. The gain (phase and amplitude) calibration was carried out in the 
normal manner after initial bandpass and delay calibration were made using the flux-density 
calibrators. All calibration and imaging steps were performed using natural weighting with 
the AIPS software package; fluxes were derived from Gaussian fits to the observed intensity 
maps, using the JMFIT task. 

Observations with CARMA (Combined Array for Research in Millimeter- wave Astron- 
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omy: a 15-element interferometer with nine 6.1 meter dishes and six 10.4 meter dishes), 
covered the three dpAGB objects AC Her, UMon & RVTau with ~ (2-3) "resolution; AC Her 
& UMon were observed at 3 and 1.3mm whereas RVTau was observed at 3mm only. Each 
source was observed in single pointings, with 15-20 min spent on source and 3-5 min spent 
on the calibrator in each observing cycle. A passband calibrator was observed in each track, 
and pointing was done every 2-3 hr for the 1-mm observations and at the beginning of each 
3-mm track. Absolute flux calibration was done using the gain and passband calibrators in 
the continuously-updated CARMA flux catalog, supplemented by observations of primary 
flux-density calibrators in those tracks in which such a source was available. Based on the 
repeatability of the quasar fluxes, we believe that our source fluxes are accurate to within ± 
20%. Calibration and imaging was done using the MIRIAD data reduction package, using 
natural weighting and the standard clean algorithm; fluxes were derived from Gaussian fits 
to the observed intensity maps, using the IMFIT task. Some of the 1.3 mm observations 
were taken in shared-risk time in which only the nine 6.1-m antennas were in the array. 



3. Results 

We detected 2 objects with the EVLA: IRAS22036 was detected in all 3 bands, whereas 
RVTau was detected in the Ka band only (Tabled]). For IRAS22036, the fluxes vary only 
weakly over the wavelength range covered by the X, Ka and Q bands, i.e., from 3.5 to 0.67cm 
(Table[T]), implying that these arise predominantly from optically-thin free- free emission. 
If, for the emitting region, we assume a typical ionized gas temperature of 10 4 K, and a 
cylindrical geometry with a length along the line-of-sight equal to its projected lateral extent, 
then the observed free-free flux, Sff, is a function of Off (source angular size), n e (the electron 
density, assumed to be uniform), and v. Thus, from the observed 33.6 GHz flux we derive 
a free-free optical depth of r//(33.6 GHz)> 0.0002 and n e > 1.1 x 10 4 cm -3 , assuming that 
Off <0"85 (i.e., the upper limit for the size of the submillimeter continuum source). By 
requiring rjy(8.5 GHz) < 1, the observed 8.5 GHz flux implies n e < 6.4 x 10 5 cm -3 and 
Off >0"067 (130 AU); the v~ 2A dependence of r ff implies 7^(33.6 GHz)< 0.055. The free- 
free emission is unlikely to arise from an ionized outflow as it would produce a spectral 
variation oc u os (Wright & Barlow 1975). 

For RVTau (not observed in X-band), our Ka-band detection and Q-band upper limit 
(107 /iJy) are inconsistent with optically-thick free-free emission and thermal dust emission, 
as these produce a spectral index, a = dlnS/dlni/ >2, implying a Q-band free-free flux 
57/(0.67 cm) ^ 490 /iJy, if scaled from the Ka-band flux. For optically-thin free-free emission, 
with a = —0.1, we expect 57/(0. 67 cm) = 262 ± 50 /iJy, consistent with our 3cr upper limit 
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of 321 /zJy. 

Mm-wave continuum emission was detected in the 3 sources observed with CARMA, 
AC Her, UMon & RVTau (Tabled]), and was unresolved (beam sizes were typically <^ 2"), 
consistent with these arising in a compact structure, likely a disk, rather than an extended 
outflow. 

If we scale the radio-to-mm/submm flux ratios of IRAS22036 for PPNe and RVTau for 
dpAGB objects, we find that the expected radio fluxes for several pAGB objects (Ka-band 
flux=2.45, 1.68, 1.54 and 1.04 mJy, respectively, for IRAS16342, 18135, 18276 and UMon; X- 
band flux=7.21mJy for IRAS19548, 20000) are far above their observed upper limits (Table 
1). This result supports the idea that the physical mechanism predominantly responsible 
for the radio emission in pAGB objects is unrelated to that responsible for the mm/submm 
emission. 

Since the observed EVLA fluxes in IRAS22036 and RVTau are insignificant compared 
to their submm/mm fluxes, the latter cannot be due to free-free emission and must be due to 
thermal emission from grains, thus providing strong support for the presence of substantial 
masses of cool large grains in these objects. Given the very low estimated free-free optical 
depth in IRAS22036, we extrapolate the observed X-band flux using a spectral index of 
—0.1, to estimate the free-free flux in the Ka and Q bands of 880 and 856 /xJy - the excess 
fluxes of 300 and 374 /xJy above the free-free extrapolated emission in these bands must 
then represent the long-wavelength tail of the large grain emission seen in the submm-mm 
wavelength range. In Fig. [I] we show the full SED of IRAS22036, together with the dust 
radiative transfer model of Sahai et al. (2006), which includes a component of cold (50 K), 
large (a = 1mm) grains. The model, which uses the DUSTY code (Ivezic et al. 1999) and 
optical constants for silicate dust provided by Ossenkopf et al. (1992), was constrained by 
the long-wavelength emission up to 2.6 mm - we find that it reproduces the excess Q-band 
flux of 374 /xJy, but the predicted Ka-band flux is a factor 3.7 times lower than our value of 
the excess, 300 /iJy. This is likely due to the inadequacy of the current dust optical constants 
for computing dust emissivity at cm wavelengths. We have therefore made new optically- 
thin dust model fits to the SED from 0.88 mm to 0.67 cm, using the Ray leigh- Jeans (R-J) 
approximation and a power-law dust emissivity, k(v) oc ^,and find that values of /5 in the 
range 1 — 1.3 provide an adequate fit ( Fig. [T] inset). 

We find that (3 = for RVTau in the 3 to 0.85mm wavelength range (computed in 
the R-J approximation), significantly lower than the corresponding value ((3 = 1.3) for 
IRAS 22036 (over the same wavelength range). If we assume that, like RVTau, the other 
two dpAGB objects, UMon and AC Her, also do not have significant free-free contributions 
at mm/submm wavelengths, then their 3 mm to 0.85 mm flux ratios yield similarly low values 
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of j3 (0 and 0.4, respectively). This difference suggests that the grains in the disks of dpAGB 
objects are likely larger than those in the dusty waist of IRAS 22036, and may be several mm 
in size: e.g., Draine (2006) finds that, for the standard power-law distribution of grain sizes 
that characterizes interstellar dust and that may apply to particles growing by agglomeration 
in protoplanetary disks, (3 <^ 1 (at A ~ 1 mm) results when grains grow to sizes a > 3 mm. 

We estimate the mass of the large grains (M^) for the 4 sources in Tabled] for which the 
observed submm flux is known to come from a compact central source, i.e, the 3 dpAGB 
sources and IRAS22036, by using the Rayleigh- Jeans approximation for grains emitting at 
wavelength A, Md = (S u X 2 D 2 )/ (2 kTdXu), where D is the source distance, and Xv is the 
dust opacity (per unit dust mass). We have not estimated dust masses for the six sources 
in Table 1 for which we are unable to associate the observed mm/submm fluxes with a 
central dusty torus (all of which are PPNe), either because the mm/submm observations 
lack adequate angular resolution (IRAS16342, 17150, 18135, 18276, 20000), or because there 
is lack of optical imaging to inform us of the object's morphological structure (IRAS19548). 

The value of Xv at submm wavelengths is uncertain by at least a factor of a few, and likely 
more. In the extensive study of the composition and radiative properties of grains by Pollack 
et al. (1994: hereafter Petal94), only large grains with a >3mm have values of < 1.3 in 
the 650/im to 2.7 mm wavelength range. If we use the largest values of Xv estimated for such 
grains from this study, then setting ^(lmm) equal to gd^t, where gd = 71.4 is the gas-to- 
dust ratio estimated from Table 2 of Petal94, and n t = 7.4 x 10~ 3 and = 1.34 is taken from 
Table 4 of Petal94, we extrapolate x„(lmm) to obtain x„(0. 85 mm) = 0.65 cm 2 g _1 . There 
are no grains in the Petal94 study with (3 as low as we have found for the dpAGB sources; 
their lowest value is (3 = 0.87 for cold 3 cm grains, for which they give x^(0- 85 mm) = 0.23 

2 —1 

cm g . 

Significantly larger values of the dust emissivity at ~ 1 mm have been used in the 
literature, e.g., Jura et al. (1997) assume x^(l-3mm) = 3 cm 2 g _1 for the carbon-rich grains 
in the PPN, the Red Rectangle. Dasyra et al. (2005) find x^(0. 85 mm) = 1.4 cm 2 g _1 from 
detailed modeling of the SED of 3 spiral galaxies. Quoting similar results from studies of (i) 
high-latitude Galactic interstellar regions by del Burgo et al. (2003), (ii) protostellar-core 
grains with ice mantles by Ossenkopf & Henning (1994), and (iii) the spiral galaxy M51 
by Meijerink et al. (2005), they argue for an enhanced emissivity at submm wavelengths. 
We derive conservative mass estimates, adopting x^(0. 85 mm) = 1.4cm 2 g _1 . The dust 
temperature is assumed to be T d = 50 K for IRAS22036 (based on the detailed modelling by 
Sahai et al. 2006), and T d = 150 K for the dpAGB sources (based on the detailed modelling 
by Gielen et al. 2007). The derived dust masses are inversely proportional to the dust 
emissivity and temperature, hence easily scaled to different values of these parameters. 
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In summary, our pilot EVLA survey confirms our hypothesis of large grain emission for 
two key pAGB objects, and thus it is likely to be applicable to these as a class - however, 
much more sensitive surveys of larger target samples at submm/mm and cm wavelengths are 
needed to establish it definitively. Current and upcoming facilities such as CARMA, ALMA 
and the EVLA will allow such surveys to be carried out with fairly modest integration times. 
We have proposed new observations with the EVLA which already allows a factor 8 larger 
continuum bandwidth (2x1 GHz compared to 2 x 128 MHz in our pilot survey), and will be 
significantly larger once the EVLA project is completed. 

RS's contribution to the research described here was carried out at the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract with NASA. Financial sup- 
port was provided by NASA through a Long Term Space Astrophysics award (to RS and 
MM). 
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Wavelength [microns] 

Fig. 1. — Observations of IRAS 22036 [black curve: ISO spectra; blue symbols: photometric 
data - squares: MSX; triangles: IRAS; circles w/center. AKARI; cross: ISO/PHOT; circles 
w/crosses: SMA 0.88 mm and OVRO 2.6 mm; diamonds: EVLA 0.67, 0.9, & 3.5 cm] and 
model spectra (magenta curve). Individual components (dashed curves) of the model are also 
shown: cool (red) & warm (cyan) shells, and hot inner disk (blue). A dashed (black) straight 
line shows a z/ -01 power-law, representing an optically-thin free-free fit to the 3.5 cm EVLA 
flux. Inset shows the long-wavelength (submm-to-cm) emission with the free-free contribu- 
tion subtracted, and optically-thin dust model fits (using the Rayleigh- Jeans approximation) 
with two values of the power-law dust-emissivity index, (3. 
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Table 1: Radio and Millimeter- Wave Fluxes of post-AGB Objects 



Source 


X 


Ka 


Q 


3 mm a 


1.3 mm b 0.85 mm 


D c 


M d 




//Jy(cr) 


yuJy(o-) 


//Jy(<r) 


mJy(cx) 


mJy(cr) mJy(a) 


kpc 


10- 2 M Q 


RVTau 




270 (50) 


(107) 


3.9(0.2) 


50.3(3.6) d 


2.2 


0.1 


UMon 




(100) 


(169) 


15(0.3) 


100(14) 182(2.6) d 


0.77 


0.064 


AC Her 


(46) 






4.6(0.4) 


38(1) 99.4(3.8) d 


1.1 


0.072 


IRAS16342-3814 


(162) 


(168) 


(254) 




277(13) c 602(90) f 






IRAS17150-3224 




(240) 


(213) 




158 (10) c 






IRAS18135-1456 


(66) 


(82) 


(169) 


12(1.4)8 








IRAS18276-1431 




(108) 


(157) 


11 (3.2)8 








IRAS19548+3035 


(45) 






6(1. l) g 








IRAS20000+3239 


(44) 






6(1. l) g 


11.4(1. 7) h 30.9(2.5)' 






IRAS22036+5306 


1010(62) 


1180(55) 


1230 (81) 


8.4(0.7) j 


290(40) j 


2 


2.2 



a Beam sizes for RVTau, UMon, & ACHcr 3mm observations are 2'/4xl'/5, 2'.'4x2'.'l, & 2'.'4xl'.'5, respectively 
b Beam sizes UMon & ACHer 1.3mm observations arc 2'.'2x0'.'9 & 2'/0xl'/8, respectively 
distances for RVTau, UMon, & ACHcr from de Ruyter et al. (2005), for IRAS22036 from Sahai et al. 
(2003) 
d de Ruyter et al. 2005 
e Giirtlcr ct al. 1996 

^adjal et al. 2010, central wavelength 0.87mm 

9 Sanchcz Contreras & Sahai 2011, central wavelength 2.6 mm 

h Buemi et al. 2007, central wavelength 1.2 mm 

l Gledhill et al. 2002 

J Sahai et al. 2006, central wavelengths 2.6 mm and 0.88 mm 



